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Abstract 

Using the QCD sum rule approach, we show that the flavor-nonsinglet H 
dibaryon states with J 71 " = 1 + , J w = 0 + , 1=1 (27plet) are nearly degenerate 
with the J 71 " = 0 + , 1=0 singlet H$ dibaryon, which has been predicted to 
be stable against strong decay, but has not been observed. Our calculation, 
which does not require an instanton correction, suggests that the Hq is slightly 
heavier than these flavor-nonsinglet Hs over a wide range of the parameter 
space. If the singlet Hq mass lies above the AA threshold (2231 MeV), then 
the strong interaction breakup to AA would produce a very broad resonance 
in the AA invariant mass spectrum which would be very difficult to observe. 
On the other hand, if these flavor-nonsinglet J=0 and 1 H dibaryons are also 
above the AA threshold, but below the S°n breakup threshold (2254 MeV), 
then because the direct, strong interaction decay to the AA channel is forbid¬ 
den, these flavor-nonsinglet states might be more amenable to experimental 
observation. The present results allow a possible reconciliation between the 
reported observation of AA hypernuclei, which argue against a stable Hq, and 
the possible existence of H dibaryons in general. 
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I Introduction 


Using the MIT bag model, JaffeJIJ predicted a stable (against strong decay) six- 
quark flavor-singlet (uuddss) hadron (referred to as H 0 ) with J=0, 1=0, and S=—2. 
A plethora of mass calculations f2| followed Jaffe’s work, and most of them predicted 
a weakly bound Hq of mass just below the AA threshold (2231 MeV). For 20 years 
experiments have searched for the H 0 , but no convincing evidence has been found 
for its existence SI- In fact, the candidate H^s from different experiments have 
very different masses. One such experiment claims the observation of a very 
weakly bound AA hypernucleus which excludes at some level the existence of the 
Ho- 

SU(3) flavor-nonsinglet H states are not usually discussed in the literature be¬ 
cause of the expectation that they should be heavier than the Ho. This expectation 
is due to the assumption that the effective magnetic one-gluon-exchange between 
the valence quarks is most attractive for the flavor-singlet channel, making the Ho 
the lightest state. 

In this work we consider the masses of two nonsinglet, doubly strange six-quark 
states which occur in the baryon-octet®baryon-octet direct product space || . These 
include the Y = 1 + , 1=0 H from the J=1 SU(3)f octet, and the 1=1, I 3 =0 H from 
the Y = 0+ 27plet. Based on model-independent assumptions, we show that it is 
likely that the H 0 is nearly degenerate with these states, and the H 0 mass is slightly 
larger. This contradiction with the MIT bag model has been discussed elsewhere 
[§]. We study the mass ratios, rather than absolute masses, because the ratio carries 
smaller uncertainty and is more stable over a wide region of the parameter space. 
The QCD sum rule approach was used by Ioffe 0] to calculate the mass splitting in 
the baryonic decuplet. 

The situation for these H states is analogous to the case of the rf (1=0, J=0), 
p (J=l, 1=1) and a; (J=l, 1=0), in which the rj' is the heaviest. The extra mass is 
attributed to the U(1)a anomalous symmetry breaking of QCD, which is taken into 
account in model calculations using instantons. The method employed here does 
not need an instanton correction because as discussed in [j§] the instanton effects (if 
any) are effectively included in the quark condensates. 

If the singlet Hq mass is greater than the AA mass, then it should not have 
been observed in previous //-search experiments ||, [|]. If the nonsinglet J=l, 1=0 
(octet) H mass is above the AA threshold (2231 MeV), but below the S°u threshold 
(2254 MeV), then it can only decay weakly or electromagnetically because the strong 
decay to AA is not allowed. The situation is different for the J=0, 1=1 (27plet) H 
since it may isospin mix with the H 0 and the J=0, 1=0 (27plet) H to form physical 
states which strong decay to AA. However, for this case, if the mixing is small, 
experimental evidence for these H s may be a narrow peak in the AA invariant mass 
spectrum. Such an observation would not contradict the observed AA hypernucleus 
events 0. The candidate nonsinglet J=0, 1=1 H s reported by Shabazian et al. 0 
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might be explained in this way. 

In Sec. P we discuss the QCD sum rule method and then formulate sum rules 
for the H 0 , the J=l, 1=0 (octet) H and the J=0, 1=1 (27plet) H. In Sec. [TIT] 
we calculate the mass ratios over a wide range of the parameter space. 

Finally, in Sec. Q we discuss our results and some experimental issues and present 
our conclusions. 


II H Dibaryon Sum Rule 


In this section, the QCD sum rules 0 are formulated for the H 0 , the J=l, 1=0 
(octet) H and the J=0, 1=1, I 3 =0 (27plet) H. We follow the method described 
in Ref. 0. The H 0 current, Jh 0 (x), is a product of two baryonic currents, Jb(x), 
which are written using the convention in [jlCj] as 


J B (x) = Of qjqk (x) + tOl' qjgk (x), 


( 1 ) 


where, 

0? mk (x) = e abc ■ [(qf a (x)Cq^(x))^ 5 q c k {x)}, (2) 

and 

Ol imk (x) = e abc ■ [(qJ a {x)C^q)(x))q c k (x)\. (3) 

In Eqs. (|2|) and (||) q is the quark held, i, j, k are flavor indices, a,b,c are color 
indices, C is the charge conjugation matrix, and the parameter t is a weighting 
coefficient for the second term of the baryonic current which is commonly used in 
QCD sum rule calculations. The Hq dibaryon singlet current is given by 

2 

Jho (x) = faijk'ti'j'k ~ -tijkti'j'k') ' JbC^$JB (4) 

where the antisymmetrization tensors, e t , jkl produce a color and flavor-singlet 6-quark 
(uuddss) state with the quantum numbers of the H 0 . 

The current correlation function is written as 


n H (q 2 ) = ~ij d 4 xe iqx (0\T J H (x) J H ( 0)|0). 


(5) 


We apply the operator product expansion (OPE) to Eq. (|5|) to obtain the right hand 
side (rhs) of the H dibaryon sum rule for large Q 2 = — q 2 . The advantage of the 
OPE for dibaryons is that there are just three non-vanishing terms that give large 
contributions in the SU(3)f limit. The resulting correlation function in the SU(3)f 

limit is 0 
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where N c is the number of color charges and (qq) is the quark condensate. The hi(t) 
coefficients in Eq. (|6|) are obtained by calculating all the possible contractions in 
Eq. (|5j). For example, for the h\{t) and t = 0, there are 3600 non-vanishing terms. 
The hi(t) s for the H 0 have been calculated elsewhere 0. 

Using the group symmetry properties of the dibaryon octet£*Doctet direct product 
space 0 the current for the J=0, 1=1, 13=0 (27plet) H is written as 

Jh 27 ,i=i,i 3 =o { x ) = Jt,o + 2J^oCq 5 J\ 

+J P Cq b J E - + J E -Cq 5 J P 

-J N Cq 5 J E o - J E oCq 5 J N , (7) 


where for the baryon currents we use the convention of Eq. (IJ). For example, 

P = ( qiqjqk ) = (udu) for the proton. The calculation yields the following hi 
coefficients: 


h x {t) = 1302 + 120t — 804f 2 + 984f 3 + 858f 4 

h 2 (t ) = 5208 + 4801 — 32167 2 + 39361 s — 64081 4 

h 3 (t ) = —20832 — 19201 + 128641 2 — 157441 3 + 256321 4 . (8) 


For the J 71 " = 1 + , 1=0 (octet) H the current is written as 


Jh*( x ) — Je+ 

JaC^Ja + \j N C^J E o + i JspC'fJN 

+J j P cy*Js- + h~-CYJ P . (9) 


This state, called the H*, does not strongly couple to the AA channel jO], |12 
this case the h, coefficients were determined to be: 


For 


hi(l) = 1866 + 20871 - 3651 2 + 14181 3 + 26361 4 

h 2 (t) = 7466 + 83491 — 14601 2 + 567U 3 — 200261 4 

h 3 {t) = —29864 — 333971 + 58401 2 — 226831 s + 801041 4 . (10) 


4 







Ill Calculation of the Mass Ratios 


The H 0 mass is given|5] by the expression 




where 

y — V S Q — sq/m 2 
* to (M*)*k\ 

accounts for the continuum part, M is the Borel mass, and So is the continuum 
threshold.^ Eq. (ill]) also holds for the other H s, so the mass ratios, given by 

r, = rn Ho (t,M^ iS^^qq) 2 ) 

m HJ (t,Mt,s { o\(qq) 2 ) ’ 

were calculated for J=0 and J=l. 

We use the standard assumptions for the phenomenological side (lhs) of Eq. (|5]) 
(pole term plus continuum for the spectral density, with continuum threshold So). 
The central values for our parameters are t = —1.2, s 0 = 5.694 GeV 2 , M— 1.5 GeV, 
(qq) 2 = (—0.250) 3 GeV 3 . We expect the calculation to be reliable for a wide range 
of the Borel mass M because the H mass is determined by the chiral symmetry 
breaking scale (~ 1 GeV)|T3| which is much larger than the scale at which QCD 
vacuum fluctuations become large (A qcd — 200 MeV) and where the Borel smearing 
fails as it does for the case of the light pseudoscalar mesons. This region is taken to 
be around 2 GeV, where the higher order terms in the OPE are strongly suppressed 
and the pole dominates the continuum contribution. Our choice for the parameter 
t, as discussed in ||, gives self-consistent QCD sum rules which do not need an 
instanton correction. Instead, the instanton effects are adequately accounted for in 
the nonperturbative vacuum condensates. The Borel mass M is in general different 
for the two H states, but we expect it to be the same for degenerate states with the 
same quark content. 

5 The third power in the Al 2 term in Eq. & is missing in Eqs. (10) and (15) of Ref. 0 due to 
a typographical error. 
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IV Results and Discussion 


Our results for the mass ratios from Eq. (0) and their sensitivities to variations 
of the input parameters are summarized in Fig. 0 . Fig- 0 a shows the dependence 
of the mass ratios on the Borel mass. The solid curve is m Ho /m H J=0 (27 P iet), while 
the dashed curve is rn,H 0 /rn,H,j =1 (octet)- The ratios remain constant, and are slightly 
greater than 1, for a large range of the Borel mass. Fig. 0b shows the dependence 
of the mass ratios on the quark condensate. The curves have the same meaning 
as in Fig. 0 a- Fig. 0b shows that large variation of the quark condensate, within a 
range consistent with that found in the literature, produces very little effect (< 2 
MeV for the mass difference). Fig. 0c and Fig. 0 d show the mass ratio dependence 
on the continuum threshold s 0 and the parameter t. Again we see little sensitivity, 
about 0.1%. Finally, the sensitivity of the mass ratio with respect to different Borel 
masses for the H 0 and the J=0, 1=1 ( 27 plet) H or the J=l, 1=0 (octet) H was 
studied. The Borel mass of the Ho was fixed and the nonsinglet H Borel mass was 
varied. These mass ratios were essentially the same at each Borel mass and they 
decreased linearly from 1.05 for a Borel mass of 1.3 GeV to 0.99 for a Borel mass of 
1.6 GeV. Other contributors to the theoretical uncertainty are the neglected terms 
in the OPE (which we expect to be small) and the vacuum saturation assumption: 

(q 2 q 2 ) ~ ( qq ) 2 , ( q A q 4 ) ~ ( qq ) 4 - 

Our conclusion, based on the results presented in Fig. |T] is that the mass ratio 
remains very close to unity in the parameter range where we believe the calculation 
to be reliable. Both the J=l, 1=0 (octet) H and the J=0, 1=1 (27plet) H are almost 
degenerate with the Ho, and this degeneracy is not sensitive to the various QCD 
sum rule parameters. The mass difference is of the order of 0.1% for both H states, 
or about 2 MeV for m#= 2 GeV. Surprisingly, we find that the singlet H 0 is slightly 
heavier than both the J=0, 1=1 (27plet) H and the J=l, 1=0 (octet) H. 

ft is worthwhile to consider the consequences of these predictions for experiments 
designed to search for neutral, strange dibaryons. We note that the nonsinglet 
J 77 = 1 + , 1=0 octet H cannot strong decay to AA (due to angular momentum and 
parity conservation), and if its mass is lower than the 5°n threshold (2254 MeV), 
then it may only decay electromagnetically (via an Ml transition to the J 71 " = 0 + , 
1=1 27plet H and/or the singlet Ho if these are lower in mass, or to the AA 7 channel) 
or weakly (if its mass is below that of the other H s and the AA threshold). If the 
J 1T = 0 + , 1=1 27plet H has mass greater than 2231 MeV it cannot strong decay to 
AA due to isospin conservation. However, isospin mixing, due to electromagnetic 
interactions within the dibaryon, cause the physical J 7r = 0 + 27plet H to contain a 
small 1=0 admixture, which allows strong decay to AA. If the admixture is small, 
then the state could have a narrow width of perhaps a few MeV. The 27plet H 
cannot electromagnetically decay directly to the singlet Ho (except via two-photon 
decay), however an El transition to the AA 7 is permitted. If the J=0 27plet H 
is below the other H s and the AA threshold, then it may only decay by the weak 
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interaction. 

If we assume the predicted mass ratios found here and that each mass is between 
the AA mass and the S°n mass, then the strong and electromagnetic decay schemes 
for the three H s will be as shown in Fig. |2j. If the decay of the J 77 = 0 + 27plet H 
is driven by the strong decay channel via the small isospin mixing 1=0 component 
rather than by the electromagnetic decay to the AAy channel, then a relatively 
narrow peak should appear in the AA invariant mass spectrum between 2231 and 
2254 MeV. The decay of the J 71 " = 1 + , 1=0 octet H to the AAy channel would result 
in a broad three-body phase space distribution for the AA mass spectrum if the 
photon is not observed. The electromagnetic decays are as follows: (1) J 71 " = 1 + , 
1=0 octet H to the J 77 = 0 + 27plet H is a AI=0 and 1, Ml transition, (2) J 77 = 1 + , 
1=0 octet H to AAy occurs via El and Ml transitions (all are AI=0), and (3) 
J 77 = 0 + 27plet H to AAy occurs via an El transition. In Fig. |]the strong decay of 
the H 0 is indicated by the thick, solid arrow, the strong decay via the small isospin 
mixing component is shown by the thick, dashed arrow, and the electromagnetic 
decays by the thin, solid arrows. 

Observation of these nonsinglet H s would not contradict the AA hypernucleus 
events already observed 0. It is possible that the candidate nonsinglet H (J=0, 
1=1) observed by Shabazian et al. 0 can be explained by these results since we 
predict the J 77 = 0 + 27plet H to be approximately degenerate with Hq. However, 
the mass must be less than that reported in |3| in order to be below the 5°n breakup 
threshold. Clearly the discovery of flavor-nonsinglet H dibaryons would require us 
to revisit traditional hadronic structure models and require a better understanding 
of quark-quark effective interactions. 

The authors thank Dr. Ron Longacre at the Brookhaven National Laboratory 
for helpful discussions. This work was supported in part by the U.S. Department of 
Energy and The Robert A. Welch Foundation. 
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Figure 1: Dependence of the mass ratio [Eq. (|T7|)1 on (a) the Borel mass, (b) the 
quark condensate, (c) the vacuum threshold, and (d) the mixing parameter. The 
solid curve is m H Jm HJ=0 (2 7p i et ), while the dashed curve is m Ho /m HyJ=1 (octe t)- 
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H Dibaryon Decay Schemes 

S N (2254 MeV) - 

H q (J=0,I=0) 

Octet (J= 1,1=0) 

27plet (J=0) 
isospin mixed 

AA (2231 MeV) 



Figure 2: Strong and electromagnetic decay schemes for the Hq, the J=l, 1=0 octet 
H , and the J=0 isospin mixed 27plet H assuming the relative masses predicted here 
and assuming the H masses lie between 2231 and 2254 MeV, as discussed in the 
text. Thick, solid arrow indicates strong decay; thick, dashed arrow indicates strong 
decay via small 1=0 isospin admixture; and thin arrows indicate electromagnetic 
decays. 
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